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Q ■ We construct an asymptotic expansion in powers of the coupling constant 

directly of the cross-section for pair production and decay of fundamental un- 
stable particles. The resonant and kinematic singularities arising in the expan- 
sion we treat in the sense of distributions. This mode allows us to transform 
formally divergent integrals into absolutely convergent ones with keeping the 
^ ', asymptotic property of the expansion. The appropriate procedure is elaborated 

' up to an arbitrary order of the expansion. The peculiarity of application of the 

procedure in the threshold region is analysed. The scheme of the calculations 
within the NNLO approximation is explicitly discussed. 

> 

^- : 1 Introduction 
o 

The processes of pair production and decay of fundamental unstable particles such 
as the top quarks or W bosons provide a powerful tool for the determination of their 
t**- ! masses and parameters of their interaction vertices. At the LHC and at an Interna- 
tional linear collider (ILC) [1, 2] the mentioned determination is planned to be carried 
out through measurings of the absolute value of the cross-section and its dependence 
^ ■ on the center-of-mass energy, and of the angular distributions and spin correlations. 

At the ILC, where the higher precision is projected, the corresponding theoretical 
calculations must be carried out in many cases with the two-loop accuracy in matrix 
elements [3, 4]. The calculation of the cross-sections thereafter must be carried out 
with the equivalent accuracy. However, the calculation of the cross-sections is ob- 
structed by resonant contributions of the unstable particles since such contributions 
are nonintegrable in the framework of conventional perturbation theory (PT). 

The common practice of a solution to this problem consists in the Dyson resumma- 
tion of the imaginary parts of self-energies in the denominators of unstable-particles 
propagators. Simultaneously this cures the Coulomb singularities, the other problem- 
atic contributions in the cross-sections [5] -[9]. However, the application of the Dyson 
resummation is rather dangerous in gauge theories because of risk to break the gauge 
cancellations necessary for unitarity. For this reason the one-loop calculations at 
LEP2 were carried out in the double-pole approximation (DPA), which guaranteed 
the gauge cancellations [10, 11]. However in the ILC case the accuracy of DPA is no 
longer sufficient [12, 13]. The application of the higher orders of the "pole expansion" 
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is ineffective, too, because of systematic inaccuracy inherent in this approach (see 
discussion in [11] and in the references therein). So in the ILC case the methods are 
required that can provide completely systematic calculations. 

Among such methods as most perspective ones long time the methods were consid- 
ered that were based on the background-field and pinch-technique formalisms. Both 
they imply a reconstruction of the conventional PT series; through the formation 
of the gauge-invariant effective action or through the Dyson resummation after ful- 
filling the gauge cancellations (see Refs. [14, 15] and [16]— [19], respectively, and the 
references therein). Unfortunately some difficulties are inherent in the mentioned 
methods, too. In particular, in the background-field formalism a dependence on the 
quantum gauge parameter remains, which cannot be fixed on physical grounds. The 
pinch-technique method does not have this problem, but there is another difficulty 
which is the common one in both cases. Namely for reaching the 0(a n ) precision 
without violating the gauge cancellation in the framework of these methods it is not 
enough to calculate corrections up to the 0(a n ) order; it is necessary to calculate all 
the 0(a n+1 ) loop corrections, too, at least their imaginary-part contributions, which 
is impractical (see, in particular, discussion in [20]). 

Another prominent approach is the "complex-mass scheme" (CMS), first consid- 
ered in [21, 22] and then worked out in [12, 13]. Actually, CMS is a renormalization 
scheme defined through the identification of the renormalized mass of an unstable 
particle with the complex pole of its full propagator. So in the CMS the renormalized 
mass absorbs the imaginary contributions of the conventional self-energy, becoming 
thus a complex quantity. Unfortunately, this triggers the complex-valued countert- 
erms, which violates unitarity. However, at the one-loop level the unitarity-violating 
contributions manifest themselves effectively at the higher orders with the gauge 
cancellations being explicitly maintained [12, 13]. So the one-loop calculations get 
legitimate. Whether the same behavior takes place at the two-loop level, this is not 
clear. For this reason the CMS cannot be considered as a rigorous procedure [13] , at 
least at present. So a development of alternative approaches still is an actual problem. 

In this paper we consider a modified perturbation theory (MPT), first proposed 
in [24] and further discussed in [25]-[28]. The basic idea of this approach is a sys- 
tematic expansion in powers of the coupling directly of the probability instead of 
the amplitude. This mode allows one to impart the sense of distributions [29] to 
the propagators squared of the unstable particles, and on this basis to asymptoti- 
cally expand the propagators squared without the appearance of the divergences in 
the cross-section. (Simultaneously, the rest of the amplitude squared is to be ex- 
panded in the conventional sense.) In doing so one gets the MPT expansion of the 
cross-section. 

It must be emphasized that the appearance of the divergences in the cross-section 
when it is considered as an integral of the expanded amplitude squared, in fact does 
not yet mean impossibility to expand the cross-section. Immediately this means only 
an invalidity of the expansion of the integrand (the amplitude squared), because 
in many cases the expansion becomes possible after the calculation of the integral. 
In such cases the distributions method is found to be a power tool to do that in 
advance, before explicit calculation of the integral. The advantage of the distributions 



2 



method is caused by its possibility to make divergent contributions convergent by 
basing on the extension principle. Unfortunately, this procedure is accompanied by 
ambiguities emerging in the integration rules. However, the ambiguities may be 
reduced to the sum composed of the (^-functions and of their derivatives with arbitrary 
coefficients, which should be added to the integrand [29, 30]. So the ambiguities can 
be eliminated via the fixing of the coefficients only. (Let us remember that the similar 
phenomenon occurs in the case of the UV renormalization at the rigid fixing of the 
renormalization scheme [31]). In the context of asymptotic expansion the coefficients 
may be fixed by means of establishing the asymptotic property of the expansion that 
should appear after the integration is carried out [30, 32, 33]. Furthermore, in some 
cases the expansion may take the form of a complete series in powers of the expansion 
parameter. At the expansion of the cross-section this would mean that each term of 
the expansion is proportional to a power of the coupling constant a and contains no 
other dependence on a. The latter property is extremely important since it means 
that in the cross-section the gauge cancellations must take place automatically at 
each power of a because the exact solution is gauge-invariant. 

At present it is generally known how to implement the MPT expansion in the 
case of pair production and decay of unstable particles, but up to the next-to-leading 
order (NLO) only and in the context of models without massless particles [28]. (This 
modelling, however, in practice has shown the existence of a smooth line-shape for the 
observable cross-section in the framework of MPT approach.) A further development 
of the method should imply a construction of the cross-section in the higher orders of 
the MPT expansion in realistic cases. Solving this problem, one should overcome dif- 
ficulties caused by the Coulomb and kinematic singularities. Recall that the Coulomb 
singularities arise as universal corrections caused by the exchanges by soft massless 
particles (photons, gluons) between outgoing massive particles in the limit of small 
relative velocities. Such contributions are nonintegrable in the cross-section in the 
higher orders of the conventional PT. However, there is a general method of summing 
up the Coulomb singularities, which makes the result integrable [6]-[9]. The kinematic 
singularities originate on account of non-analyticity of the phase-volume factor. The 
non-analyticity, in turn, gives rise to nonintegrable singularities in the cross-section 
in the presence of the MPT-expanded propagators squared. Consequently the latter 
singularities is an essential problem in the MPT approach. 

In this paper we find a solution to the problem of kinematic singularities in an 
arbitrary order of the MPT expansion. The Coulomb singularities are taken into 
consideration rather qualitatively, in the spirit of [6]-[9], which is sufficient for our 
purposes. We emphasize that we search for a principle solution to the problem of 
curing non-integrable contributions in the cross-section. Correspondingly, we carry 
out only analytical calculations necessary for reducing further calculations to the 
realizable ones by means of numerical methods. In doing so we mean that the con- 
ventional regular contributions to the cross-section can be independently calculated 
in the framework of conventional PT. 

The paper is organized as follows. In Sect. 2 we formulate the problem and do 
its preparatory analysis. An algorithm of the calculation of the MPT expansion 
in the case of the hard-scattering cross-section, is worked out in Sect. 3. In Sect. 4 
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we calculate singular contributions to the observable cross-section and introduce a 
necessary modification to the method in the near-threshold regions. In Sect. 5 we 
discuss the peculiarities of the calculations within the NNLO approximation. Sect. 6 
outlines the results. In Appendix A we derive a formula for the characteristic change 
of variables in the (^-function with derivatives. In Appendix B we carry out analytical 
calculations of basic singular integrals. 

2 Statement of the problem 

We analyse pair production and subsequent decay of unstable particles, in particular 
in the process of e + e" annihilation. For simplicity we consider the case when only 
stable particles are immediately produced at the decays of the unstable particles. (The 
cascade decays should not lead to new significant difficulties.) The observable cross- 
section of the whole process is determined via the convolution of the hard-scattering 
cross-section with the flux function describing contributions of nonregistered photons 
emitted in the initial state [10]: 

s 

f ds' 

a(s)= J —4 > {s'/s-s)a{s'). (2.1) 

Here s is the energy squared in the center-of-mass system, a(s) and a(s) are the 
observable and the hard-scattering cross-sections, respectively. Apart from s both 
cross-sections can depend on spin and angular variables. (We do not consider explic- 
itly this option as the corresponding modifications are unessential for our analysis.) 
The ratio z = s'/s characterizes a fraction of the energy expended on the production 
of unstable particles. For our purposes it is sufficient to take flux function <f>(z; s) in 
the leading logarithm approximation. So we put 

<f>(z; s) = (3 e (l - z)^-V - I/3 e (l + z ), (3 e = ^ L -L _ lY ( 2 .2) 

2 it \ m 2 e J 

The hard-scattering cross-section a(s) must be determined in the PT framework, 
though in some cases indirectly because of divergences arising in the phase-space 
integral. To avoid the divergences, usually the Dyson resummation is applied in 
the unstable-particles propagators. Starting conceptually from an exact solution, we 
initially imply the application of the Dyson-resummation, too. However further we 
will asymptotically expand the propagators squared in the sense of distributions. 

Generally, there are two types of contributions to the hard-scattering cross-section, 
the factorizable and non-factorizable ones. The factorizable contributions retain the 
structure of the process as the sequential production and decay of unstable particles. 
The non-factorizable contributions connect these subprocesses. Typically the non- 
factorizable contributions are represented by configurations with contributions of one 
of the unstable particles or of both of them being components of one-particle irre- 
ducible subdiagrams. Such contributions do not generate singularities in the ampli- 
tude and therefore are integrable in the cross- sect ion. Generally, the non-factorizable 
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contributions can be related to the cases of inclusive single production and decay 
of unstable particles or direct production of final states. The corresponding calcu- 
lations in the MPT framework are simpler as compared to the cases with factoriz- 
able contributions. For this reason we postpone the consideration of non-factorizable 
contributions to the end of the next section. The only exception we make is for 
soft massless-particles contributions. The matter is that they can retain the double- 
resonant structure even if they formally are non-factorizable [34]-[36]. (Actually this 
property is a manifestation of general theorem [37, 38] which states that the soft 
massless-particles contributions are collected in a factor in the cross-section.) For 
this reason we will consider the soft massless-particles contributions in parallel with 
the consideration of factorizable contributions. 

So, at first we consider the piece of the cross-section that possesses the double- 
resonant structure. The hard-scattering cross-section in this case is conveniently 
written as follows: 

oo oo 

o-(s) = JJds 1 ds 2 9(\ / s-y/s^-^)a(s;s 1 ,s 2 )(l+5 c ) (2.3) 

5 lmin 
5 2min 

Here a(s ; s±, s 2 ) is an exclusive cross-section, the (1 + 8 C ) stands for contributions of 
soft massless particles, Si and s 2 are the virtualities of the unstable particles. The 
si m i n and S2min are the minimums of Si and s 2 , respectively, defined as squared sums 
of the masses of the decay products of the unstable particles. Hereinafter we imply 
that the square roots of si m i n and s 2m in are less than the masses of the corresponding 
unstable particles. The exclusive cross-section a(s ; s±, s 2 ), we write in the form with 
extracted Breit-Wigner (BW) factors: 

a(s ; s h s 2 ) = \y/X(s,Si,s 2 ) $(s; s u s 2 ) pi(si) p 2 {s 2 ) . (2.4) 

Here Pi(si) are the BW factors, A(s, s±, s 2 ) = [s—(^/sl+^/s 2 ) 2 ][s—(^/sl—y/s 2 ) 2 } is the 
so-called kinematic function. The product of the a/ A(s, s±, s 2 ) with the ^-function 
in (2.3) constitutes the kinematic factor. Let us note at once that the kinematic 
factor has a singularity at yHs = y/s\ + y/s^. The function $(s; s\, s 2 ) is the rest of 
the amplitude squared. By construction it does not include singularities associated 
with the production and decay of unstable particles. Correspondingly, we consider 
$ determined in the framework of conventional PT. The BW factors, we define as 
follows: 

P^H^xIA^)! 2 . (2.5) 

7T 

Here Mj is the renormalized mass and Tj is the Born width of the ith unstable particle, 
Aj is a scalar part of the Dyson-resummed propagator (so the spin factors are referred 
to $). Further we omit index V if its presence is not necessary. The A(s) has the 
form 

A(S) = s - M 2 + ReS(s) + i ImS(s) ' (2 ' 6) 
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where ReS(s) and ImS(s) are the real and imaginary parts of the renormalized self- 
energy. We assume that S(s) is determined in the on-mass-shell (OMS) scheme of the 
UV renormalization, or rather in some version of its generalization to the unstable- 
particle case [39]-[43]. Recall that the naive expansion of A(s) in powers of the 
coupling leads to divergences in integral (2.3). 

Now we proceed directly to arguing the MPT approach. Its most important feature 
is the assignment of the distribution sense to the factors that explicitly or potentially 
include singularities. In the given case such factors are the BW factors, the kinematic 
factor, and the factor 1 + 8 C . Function $ may be considered as a representative of 
trial functions since $ is smooth in the vicinity of prospective singularities (located 
actually on the thresholds). By the next step, it is necessary to asymptotically expand 
the product of the mentioned factors in the sense of distributions. 

For methodological reasons at first we consider the asymptotic expansion of the 
B W factors by considering each of them isolated. A general structure of their expan- 
sion is as follows [24]: 



N 



p(s) = 5(s-M 2 ) + PVT N { p(s) } + c «(«) $n(s-M 2 ) + 0{a N+1 ) . (2.7) 

n = 

Here the expansion is carried out in powers of a up to neglected terms of order 
0(a N+1 ). The leading term originates by virtue of the relation a/ (x 2 + a 2 ) — > nS(x) 
as a — > and the unitarity relation ImS(M 2 )| 1 _i oop = Mr|B m- The second term 
denotes the Taylor (T) expansion in powers of a up to with the poles in s — M 2 
defined in the sense of principal value (PV). Recall that the principal value may be 
defined by the relation 

1 d n 
PV x™ ~ (ra-1)! dx™ ln ^ ' 

where the derivatives are understood in the distributions sense, i.e. at an integration 
they should be switched to the weight function via formal integration by parts. The 
sum in (2.7) corrects the contributions of the PV-poles in the points where the PV- 
prescription operates. The S n (. . .) denotes the nth derivative of the 5-function taken 
with a typical coefficient, 

(—\ n 

5 n (x)= K -^-5^(x). (2.9) 

The highest degree of the derivatives of the 5-functions in (2.7) is equal to the highest 
degree minus 1 of the PV-po\es in the Taylor expansion. The coefficients c n (a) 
are defined so that to provide the asymptotic property for the expansion. In this 
way the contributions of ^-functions in (2.7) mean finite counterterms. As for c n (a) 
themselves, they are the polynomials in a with the exponents ranging from n to N. 
The coefficients of the polynomials are defined by self-energy of the unstable particle 
and by its derivatives calculated on the mass shell. (The contributions of soft massless 
particles are considered regularized at this stage.) 

For foreseeable applications, seemingly, the expansion within the NNLO approxi- 
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mation will be sufficient. With this precision formula (2.7) takes the form 



9N MT f 1 2aReS!(s)l 

, (S ) = ,(._«>) + [PVj^j^ - PV 1I - J AI\ (2.10) 

2 

+ ^c n (a)5 n (s-M 2 )+0(a 3 ). 

ra = 

Here Si is the one-loop self-energy. The /c-loop self-energy S fc is defined through the 
formula 

= a Ei(s) + a 2 S 2 (s) + a 3 E 3 (s) + • • • . (2.11) 

The polynomials c n (a) within the NNLO approximation actually have been calculated 
in [24], though have been presented not quite explicitly. Below we represent c n (a) 
in completely expanded form in an arbitrary version of the above mentioned gener- 
alization of the OMS scheme. All these versions are characterized by the one-loop 
conditions Re£i(M 2 ) = 0, Re£^(M 2 ) = 0. So we have 

co = -«| + « 2 (|-f + (2.12) 

d = - a 2 (hl( + R 2 ) , c 2 = -a 2 I 2 . 

Here I k = Im£ fc (M 2 ), R k = ReS fc (M 2 ), and the primes mean the derivatives at 
s = M 2 . The i?2 and R' 2 are determined by the renormalization conditions, too, 
though differently in different versions of the generalization of the OMS scheme [42]. 
In particular, in the OMS or "pole" scheme [42, 43], R 2 = -IJ{ and R' 2 = -IJf/2. 

So, we have determined the asymptotic expansion of isolated BW factors. Now 
we consider the product of BW factors. It is clear that in the case of smooth weight 
the asymptotic expansion of the product has the form of the formal product of the 
expansions of separately taken BW factors. However, this is not the case if the 
weight includes singularities which intersect with the singularities of the expanded 
BW factors. Unfortunately, this occurs in the case of integral (2.3) because of the 
soft massless-particles contributions and the singularity of the kinematic factor. In 
order to asymptotically expand, nevertheless, the integrand in (2.3), we should in- 
dependently make definition of the expansion of the product of BW factors in the 
vicinities of intersections of singularities. 

At first we discuss the case of the singularities caused by the soft massless-particles 
contributions. Recall that we symbolically expressed them through the factor 1 + 5 C 
in formula (2.3). Actually these contributions become crucial when the soft massless- 
particles momenta merge with the momenta of the unstable particles (owing to the 
emission or absorbtion of massless particles) and when the unstable-particles momenta 
turn out to be on the mass shell. The asymptotic expansion in this case can be fulfilled 
if the additional singularities, intersecting with the singularities of the BW factor, are 
cancelled in the cross-section. In reality the cancellation takes place in the case of 
the ordinary IR-divergent contributions, but not in the case of Coulomb singularities. 
Therefore a solution to the problem of soft massless-particles contributions can exist 
in the case only when Coulomb singularities are somehow regularized. 
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Let us postpone the problem of the Coulomb singularities and concentrate on the 
ordinary IR-divergent contributions. A solution to this problem actually is proposed 
in [24]. It is based on the introduction of dimension regularization for the IR sin- 
gularities and then on the calculation of additional counterterms at the intersections 
of singularities. Apart from the (^-functions of virtualities of the unstable particles, 
these additional counterterms include also the (^-functions of the massless-particles 
momenta. Moreover, the coefficients at these counterterms are found to be singular 
in the dimension-regularization parameter. However the later singularities are to be 
cancelled after the integration over the massless-particles momenta. 1 In this way the 
counterterms become identical in form with those in the case without the contribu- 
tions of soft massless particles (though they appear with modified coefficients). An 
alternative solution is based on the introduction of the IR regularization through 
inserting of a technical mass for massless particles. A remarkable property of this 
regularization is the absence in its framework of the above mentioned intersections 
of singularities [26]. By virtue of this fact the extra counterterms do not arise at 
all in this case. So, again, after the integration over the massless-particles momenta 
the structure of the expansion becomes the same as in the case without the contri- 
butions of massless particles. For this reason we further do not consider explicitly 
the contributions of soft massless particles except contributions that lead to Coulomb 
singularities. 

Proceeding to the Coulomb singularities, we note first of all that they arise in a 
more restricted kinematic region and have different nature as compared to the or- 
dinary IR-divergent contributions. So from the very beginning they can be isolated 
from the ordinary soft massless-particles contributions and thereupon can be sepa- 
rately investigated. By this means it was found that the Coulomb singularities arise 
on the mass shell and simultaneously in the limit of small relative velocities of out- 
going massive particles. At the level of the cross-section they may be collected in the 
so-called Coulomb factor. For the first time this factor was discovered in the case of 
stable particles and in electrodynamics with the result [44, 45] (hereinafter we omit 
the ordinary IR-divergent contributions) 



Here ^(0) is the wave function at the origin of the charged massive particles moving 
in the center-of-mass frame (c.m.f.) with velocity /3q, k is a positive coefficient. (In 
the case with color particles in QCD the Coulomb factor is the same with the modi- 
fication in the coefficient k only [46, 47]). In the expanded form factor (2.13) can be 
represented with the aid of the formula 



1 The basic part of the calculation of the additional counterterms in the case of pair production 
and decay of unstable particles, at the level of the one-photon contributions, has been made by the 
author of this article in collaboration with F.Tkachov [25]. 



1 + 5 C = |^(0)| 2 



X 




(2.13) 





(2.14) 
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where B n are Bernoulli numbers. 

Returning to the MPT expansion, we note that the PV-poles and the derivatives of 
the (^-functions in the expansions of the BW factors are ill-determined in the presence 
of the Coulomb singularities. Really, in accordance with the integration rules of the 
PV-poles and the derivatives of the (^-functions the Coulomb factor must be taken 
into consideration not only on-shell, but also with the derivatives with respect to 
virtualities of unstable particles. However the derivatives of the Coulomb singularities 
cannot be determined as they appear precisely on the mass shell. 

In this connection further we use an extended representation of the Coulomb factor 
available in the unstable-particle case. Namely we mean a particular resummation of 
the singular and some kind of nonsingular Coulomb contributions, which should be 
made at a given number of exchanges by soft massless particles. The resummation 
should be made for the off-shell outgoing unstable particles and, simultaneously, with 
the Dyson resummation in those propagators of the unstable particles that are directly 
involved in the generation of the Coulomb singularities. 2 The practical extraction of 
the mentioned contributions is a rather subtle task. For details we refer to [6]-[9], 
where the extraction is made both for one- and multi-photon contributions. Below 
we adduce the result with the one-photon contribution. In the notation of [8], but 
omitting some superfluous term [9] which is inessential for the validity of the ultimate 
result, we have 

na 

c | one-photon — 2f3 

si, s 2 ) is the velocity of outgoing unstable particles in the center- 
of-mass frame, (3m = a/1 — 4(M 2 — iMT)/s, T is the width of unstable particles. 
Notice that with nonzero T the r.h.s. in (2.15) is a smooth function of (3. At the 
same time at going to the mass-shell and subsequently taking the limit r — > 0, the 
r.h.s. becomes the Ka/(2(3 ), which is the one-photon Coulomb singularity. 

We emphasize that in fact only imaginary part of the on-shell self-energy may be 
involved in the above resummation. Therefore the regularized Coulomb factor may 
be made gauge-invariant in the OMS or "pole" scheme of the UV renormalization 
[42, 43], because in this scheme both the renormalized mass and the imaginary part 
of the on-shell self-energy are gauge invariant quantities. The gauge invariance of the 
regularized Coulomb factor, in turn, implies the gauge invariance of the piece of the 
cross-section that remains after the extraction of the Coulomb factor. So provided 
this remaining piece is completely expanded in powers of a, the contributions at each 
power of a must possess the property of gauge cancellations. In the next section we 
show the existence of the complete expansion of the above mentioned remaining piece 
of the cross-section. 

So we have gained the regularity property of the Coulomb factor and the possi- 
bility to differentiate it with respect to the virtualities of the unstable particles. In 
addition, we have a regularity of the derivatives of the Coulomb factor, too. The lat- 
ter property most easily may be shown by noting that the regularized Coulomb factor 

2 The mentioned propagators belong to one-particle irreducible subdiagrams. Therefore they by 
no means are involved in the BW factors that are to be MPT-expanded. 
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(2.15) 
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is an even function of (3. In the multi-photon case this is evident from the appropriate 
explicit expression [9]. In the one-photon case this becomes clear if proceeding to the 
representation 

*c| on e-photon = 5 ^(l^F^) ' ^ 
which is equivalent to (2.15) at small enough /3. From the even property it follows that 
S c depends effectively on f3 2 . (In particular this is evident from the Taylor expansion.) 
This implies that both the 5 C and the derivatives of 5 C with respect to are smooth 
functions of Sj. 

The property of smoothness of the regularized Coulomb factor and also of its 
derivatives with respect to the virtualities allows us to refer the Coulomb factor to 
the weight in formula (2.3). Correspondingly, at the further calculations we consider 
factor 1 + 5 C as absorbed by function $. In doing so, the 1 + 8 C may be taken 
into consideration in the expanded form with respect to the number of exchanges by 
soft massless particles. Unfortunately, this is not the complete expansion since some 
dependence on a is involved in the parameter of the regularization. Nevertheless with 
the fixed parameter T, the expansion is an asymptotic one with respect to a. 

Finally, we turn to the kinematic factor 9(y/s — ^fs\— )a/A(s, ~s\, 82). As noted 
above, its singularity located on the point set v^i — v^2 = can intersect with the 
singularities of the expanded BW factors. Most easily this may be shown by imposing 
conditions on s. So, at y/s = v /sJ + M; L (i / j) the point set y/s — yfsl — y/s2 = is 
reduced to the isolated point yfsl = Mi, where the singularity of one of the expanded 
BW factors is located. At y/s = M1 + M2 the above point set is reduced to y/s~i+y/s~2 = 
Mi + M2, where the singularities of both BW factors are simultaneously located. 

At first sight this property makes senseless the rules of the integration of the 
PV-poles and of the derivatives of the 5-functions, and hence the very possibility of 
the calculation of the asymptotic expansion of integral (2.3) becomes questionable. 
However, the distribution method in this case again allows us to add a sense to 
the integrals that arise in the course of the expansion. Moreover, the ambiguities 
which are peculiar to this method may be again removed by means of the imposing of 
the asymptotic condition of the expansion. Furthermore, both these requirements are 
automatically satisfied by application of the analytical regularization of the kinematic 
factor. In the next section we proceed to a systematic discussion of this issue. 



3 The scheme of calculation of &(s) 

Fist of all we reformulate the problem in dimensionless variables. For this purpose 
we introduce dimensionless energy variables x, xi, X2 counted off from thresholds, 

V^=2M + yX, y/s~= Mi+^-Xi. (3.1) 

Here i = 1,2, and M = (M 1 + M 2 )/2. Then formula (2.3) takes the form 

a(x) = / / dxi dx 2 9{xi + ai)9{x2 + a2)9{x—xi — X2)\/x — xi—X2 (3.2) 
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$(x;xi,x 2 ) [l + 5 c ) pi(x 1 )p 2 (x 2 ) 



Here a(x) = M 1 M 2 a(s), Pi(xi) = MMiPi(si), and a* = 2(Mj - y/s^~)/M. At once 
we notice that the parameters a* are strictly positive at > Sj min and that in any 
case ai + a 2 < 4. Notice also that in view of the ^-functions, the a(x) is effectively 
proportional to 9(x+ai+a 2 ). The 1 + 5 C in (3.2) is equivalent to within the change of 
variables to 1 + 5 C in (2.3). Functions $ and $ are related each other by the formula: 

y 1 

dxi dx 2 y/x-x 1 -x 2 $(x ; x 1 ,x 2 ) = 777 dsi ds 2 -7 a/A(s; si, s 2 ) $(s; si, s 2 ) . (3.3) 

Since the factors y^ - x i — ^2 and a/ A(s; s\, s 2 ) have identical analytical properties 
in the integration area, the analytical properties of $ and $ are identical in this 
area, too. Let us remember that by the construction $ is smooth in the vicinity of 
physical thresholds associated with the production and decay of unstable particles. 
So $(s; si, s 2 ) has no singularities in the points y/s = 2M, yfsl = Mi, \/s — yjs\ m \ n — 

\J s 2 min = 0, — yj Si m ; n — Mj = 0, y/s — yj Si m ; n — ^/~Sj = 0, — yJ~S~i — Mj = 0, 

\/s — V^i — V^ 2 " = (i 7^ j). Correspondingly, $(x ; xi, x 2 ) has no singularities in the 
points x = 0, Xi = 0, x + a = 0, x + a>i = 0, x + a>i — Xj = 0, x — xt = 0, x— X\— x 2 = 0. 

The asymptotic expansion of isolated BW factors Pi(xi) most easily may be ob- 
tained by means of change of variables in formula (2.7). So, let us rewrite the second 
relation in (3.1) in the form s t - Mf = M.\I,.r, \ + XjM/(4Mj)], and note that 
[1 + Xi M/(4Mj)] ^ everywhere in the integration area. This implies that the sin- 
gular point Si — Mf = transforms into the point Xi = in new variables. So the 
5(si~ Mf) becomes the (MM i )~ 1 5(x i ), while the derivatives of S(si-Mf) transforms 
into the sum of S(xi) and its derivatives. The final formula valid within the integration 
area is derived in Appendix A, 



Recall that the 5 n means the ra-th derivative of the 5-function with a typical coefficient, 
cf. (2.9). As long as the PV-regularization is a canonical one [30], the change of 
variables in the PV-poles is made as in conventional functions, i.e. without the adding 
of a sum of the 5-function and its derivatives: 



Substituting (3.4) and (3.5) into (2.7), we get the asymptotic expansion of Pi(xi). 

Unfortunately, the presence of the kinematic factor 6{x — x\ — x 2 )y/x — x\ — x 2 
does not allow us to substitute the expansions for p~i(xi) in formula (3.2), because 
the derivatives of the kinematic factor are non-integrable, or even completely un- 
determined in the point {x — Xi = 0, Xj = 0}, i ^ j. Nevertheless the situation 
changes if to proceed to the analytically regularized kinematic factor, namely to 






(3.5) 
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6(x— x\— x 2 ) (x—x\—x 2 ) x with non-integer A. 3 With large enough A this gives a regu- 
larity property for the derivatives. As a result, after the substitution of the expansions 
for Pi(xi), the integral (3.2) considered as iterated becomes calculable. Further, after 
the integral is calculated, we have to take the limit A — > 1/2 in order to match the 
original integral. So the problem is reduced to the taking of the limit of the already 
calculated integral. Notice that the above procedure means the treatment of the 
kinematic factor as the distribution (x — X\ — x^ 1 / 2 [30]. 

Fortunately, if the above mentioned limit exists then this guarantees the asymp- 
totic property of the expansion of the whole of integral (3.2). This follows, first of 
all, from the fact that prior to any expansion the substitution of the distribution 

1 /2 

(x— X\— x 2 ) + for the kinematic factor does not mean a modification of the integral. 
This, in particular, means that the asymptotic expansion after the substitution is a 
completely rightful operation. Technically the expansion may be carried out through 
the proceeding at intermediate stages of calculations to the large enough A. (The 
infimum of A is dependent on the order up to which the expansion is to be made.) 
This allows us to treat the regularized kinematic factor as a smooth weight and on 
this basis to asymptotically expand the BW factors. The integrals that will appear 
after the expansion by our construction will be calculable. Further, the analyticity 
will allow us to vary A with keeping the asymptotic property of the expansion. The 
existence of the limit at each term of the expansion will imply the implementation of 
asymptotic property of the expansion. 

Now let us proceed to the calculation of the above mentioned integrals at arbitrary 
but large enough A. For this purpose, we consider any one term of the resulting 
expansion of the integrand in (3.2) that appears after the formal expansion of the BW 
factors. Let this term include PVx± ni or 5 ni -i(xi) from one of the BW factors and 
PVx 2 n2 or 5 n2 -i(x 2 ) from another BW factor (711,2 > 1)- Recall that we consider the 
kinematic factor in the form (x— x\— £2)+ having in mind that at the intermediate stage 
the A is large enough. The Q{x\ + ai) 9(x 2 + a 2 ), we consider as a potentially singular 
factor, as these ^-functions via the integration by parts can generate non-integrable 
singularities at A = 1/2. As a pure weight we consider the function $(2; ; x\, x 2 ) 
multiplied by the regularized Coulomb factor and by the definite regular coefficients 
that arise in the course of application of the formulas (2.7), (3.4) and (3.5). For 
brevity we denote the mentioned weight by the same symbol $(x ; xi, x 2 ). 

Further, let us subtract and add to $ the n\ x n 2 of the first terms of its Taylor 
expansion in powers of x 1 and x 2 , 

ni-ln 2 -l ki k 2 _ 

=EEnF! $(tl,h) ( i;0 ' )- (3 - 6) 

fc 1= 0fc 2 =0 1- 2 ' 

Then $ takes the form of a sum of two terms, the first one is (3.6) and the second 

3 Let us remark that (x— x\— X2) x instead of y/x— x\— X2 arises automatically if the phase- volume 
integration is carried out at the space dimension d = 4 + 2e,e = A— 1/2. 
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one is the difference A$ = $ — The difference, we represent in the form 

"i- 1 r ki "2-1 k 2 

a$(x ; xi, x 2 ) = E + E n 

fei=0 1- fc 2 =0 2 ' 

+ i; i ^>(x;s 1 ,i 2 ). (3.7) 

Here the first sum represents the Taylor expansion up to of the remainder of 

the Taylor expansion of $(x ; x i: x 2 ) in x 2 up to x^ 2 ^ 1 . The second sum has the 
same meaning with the replacement x x <-> rr 2 . The third term has the meaning of 
the common remainder. Note that by the construction the functions ipn*\x,Xj) and 
ip(x ; xi, x 2 ) are regular in some neighborhoods of the points x\ = 0, x 2 = 0. 

It is clear that leads to singular contributions in the presence of the PV- 

poles and the (^-functions. However, thanks to its power-like dependence on x\ and 
x 2 the integral (3.2) of these contributions can be analytically calculated. Namely, 
owing to the relations 

(y-x) x + x k 5 n (x) = (y-x) x + 5 n _ k (x), 0<k<n, (3.8) 

{y-x)\x k PV^ = (y-x) x + PV^:, 0<k<n, (3.9) 

which are valid in the framework of the analytical regularization, the mentioned 
contributions to (3.2) are reduced to linear combinations of the basic integrals 

A k- L k 2 ( x ) = JJdx!(lx 2 9(x 1 + a 1 )9(x 2 + a 2 ) (x-x 1 -x 2 )+5 kl - 1 (x 1 ) S k2 -i(x 2 ) , (3.10) 
B x ik2 {x) = JJdx 1 dx 2 9(x 1 +a 1 )9(x 2 +a 2 )(x-x 1 -x 2 ) x _PV^S krl (x 2 ), (3.11) 

C k \ k2 {x) = I Ux 1 dx 2 9{x 1 +a 1 )9{x 2 +a 2 ){x-x 1 -x 2 ) x + Pv\- i Pv\. (3.12) 

Here subscripts fcj range over 1 < ki < rij. Notice that among the S-type inte- 
grals that include both the PV^-pole and the 5-function, there are integrals with the 
opposite order of PV and 5. Let us denote these integrals by B k x ik2 . We do not 
write down them explicitly as they may be obtained immediately by substitutions 
{ki,ai} <-> {^2,02} in formula (3.11). 

The analytical calculation of integrals (3.10)-(3.12) is carried out in Appendix B. 
In the general case the outcomes have the form of a sum of regular and singular 
at A = 1/2 contributions. The singular contributions are defined as the power-like 
distributions of the type x\ multiplied by some regular factors, see Table 1. The 
calculation of the observable cross-section in the presence of the singular distributions 
is considered in the next section. Here we remark only that with a complicated 
weight the integral of the distribution x+ cannot be analytically calculated, but can 
be calculated numerically with the aid of formula [30] 

ldxx x + ^(x)= I dxx x lip(x) - |[^ (fc) (°) f • ( 3 - 13 ) 
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Table 1: The singular contributions to a(x) arising via the two-fold basic integrals (r 
is a no n- negative integer such that 5/2— n 2 + r < 0). 



A 1/2 


5/2-ni-m 
X + 




Rl/2 




/ . \5/2— n2+r 

(x + ai) + ' 




5/2— ni— r«2 

x + 





Here A < — 1 and A is non-integer, K is a positive integer such that — K— 1 < A < —if, 
?/?(x) is a test function which must be smooth enough in the vicinity of x — 0. 

In the difference term A<E>(x ; Xi, x 2 ), we consider separately the contributions of 
the sums and of the common remainder. In the common remainder in view of the 
relations 

x n 8 n - 1 {x) = 0, x n PVx~ n = l, (3.14) 

the multipliers x™ 1 and x 2 2 make vanish contributions of the 5-functions and cancel 
the PV-poles. In order to calculate the remaining integral, we proceed to the cone 
variables x\ + x 2 = £, X\ — x 2 = 2r] and make a shift £ — > rr — £. After the calculation 
of the integral d.77, the relevant contribution to (3.2) takes the form of 

J d£ 0(x + ai + a 2 - £ ^(x, . (3.15) 

Here (p(x,£) is a regular function that arises after the calculation of the integral drj. 
In the case A = 1/2 integral (3.15) is absolutely convergent and therefore can be 
numerically calculated. 

In the sums in A<E>(x ; Xi, X2), the multiplier x" J makes vanish contributions of 
8nj-i(%j) an d cancels PVx~ Uj . The contributions of 5 nj _i(x i ) and PVx~ ni remain 
not touched upon, but the dependence on Xj is simpler. Therefore at first we calculate 
the integral dxj. In view of (3.8) and (3.9), it is reduced to the sum of the basic 
integrals (1 < k < rij) 

I^(x - xj) = j dxi9(xi + di)(x - Xi- Xj)\8 k -i{xi) , (3.16) 

J x k (x - Xj ) = / dx, d( Xl + ai ) (x-Xi- Xj)\ Pv\ . (3.17) 
J %i 

The analytical calculation of these integrals is carried out in Appendix B, as well. 
Having analytical expressions for \\ and j£, we can calculate the remaining integral 
dxj . In the general case it is reduced to the sum of integrals of the kind 

J dx j 9(x j + a j )^;\x,Xj)P x k (x-Xj). (3.18) 
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Table 2: The singular contributions to a(x) associated with the one-fold basic inte- 
grals. (Here r is a non- negative integer such that 5/2 — k + r < 0.) 

4 12 (-* - <*)P -' + ' 



Here stands for l£ or j£ and ifnj\x, Xj) is introduced in (3.7). By making the 
change of variable x — Xj — > £, we rewrite (3.18) in a more convenient form, 



d£ 9(x + a, - ^nT(x, x - P£(0 • (3.19) 

At first we consider the case = l£. In view of l£(0 ~ £+ +A_/ \ integral (3.19) 
is absolutely convergent at A > A; — 2. At A < A; — 2 and x + Oj ^ the integral can 
be calculated with the aid of (3.13). The condition x + aj ^ is necessary for the 
providing of smoothness for the weight at £ = 0. As x + aj — > 0, integral (3.19) tends 
to infinity. Nevertheless, the singular contribution can be analytically calculated and, 
simultaneously, the regular background can be numerically calculated. Namely we 
represent ifn-(x,x — £) in the form of the Taylor expansion with a remainder, 



K-1 

<pSt\x,x- = E 7[ rf5" r W) + A^>(x,x - • (3.20) 



r=0 



Here -X-1 < 1 + A-A; <-K, and (p^' r \x, x) = d r /dC (p%?\x, x-£)\ e=0 - Owing to 
Atpnj\x, x — £) = 0(£ x ), the integral with A(pn^ is absolutely convergent, and this 
contribution can be numerically calculated. The integral of the sum in (3.20) gives 



E 



2 + A-A; + r r! ' 1 ' ' 

r=0 

Here 2+A— A;+r < 0, and we note the distribution sense of the singular contributions. 
Up to a coefficient, formula (3.21) determines the appropriate singular contributions 
to a(x). 

In the case = j£, we have generally two singularities in P£(£), the (— £)+ +A ~ fc 
and £+ +A_fc , cf. (B14). However, the latter singularity does not make contributions 
at A = 1/2 in view of the factor cos(7rA). To calculate the contribution of the former 
singularity, we do the change of variable £ — > — £ and take advantage of the relation 
9(x + cij +£) = 1 — 9{— x — aj — £). It is clear that the unity gives regular contribution. 
The 0(— x — aj— £) gives a regular contribution outside of a neighborhood of x+aj = 0. 
Inside the neighborhood, by complete analogy with the previous case, it gives a sum of 
the regular contribution and the power- like singular contributions (— x — a ? -)^ +A ~ A:+r . 
The singularities arising in both cases are shown in Table 2. 
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Thus, we have reduced integral (3.2) to the sum of regular and singular contribu- 
tions. The regular contributions are described by the absolutely convergent integrals 
and thereupon can be numerically calculated. The singular contributions are repre- 
sented by the power-like distributions with regular weights. So their further integra- 
tion can be carried out with the aid of formula (3.13). It is worth noticing that the 
singularities that remain at the level of hard-scattering cross-section a(x) are located 
on the physical thresholds x = and x = — a^. 

So, we have discussed the piece of the hard-scattering cross-section that initially 
possessed the double-resonant structure. However, this piece includes not only the 
double-resonant contributions but through the difference term (3.7) also the single- 
resonant and the non-resonant contributions. Now we consider the piece of the cross- 
section which initially possesses the single-resonant structure and which is originated 
by the integral of the type (3.2) without one of the BW factors. The kinematic factor 
in this case is present as before, but one of its variables, say, Xj describes immediately 
the invariant mass of the appropriate particles in the final state. On the contrary, 
the Xi variable describes first of all the virtuality of the unstable particle and only 
then the invariant mass of the appropriate stable particles. In this case we do at first 
the MPT expansion of the BW factor which depends on Xi variable. Then the weight 
at each term of the resulting expansion, we represent in the form of a finite power 
series in x^ with a remainder giving regular contribution. The latter contribution may 
be reduced to an integral of the type (3.15), which is to be numerically calculated. 
The integral dxi of each term of the finite power series, we reduce to the sum of 
the basic integrals (3.16) and (3.17). After that the integrals dxj we reduce to the 
sums of the integrals of the type (3.18). So we get again the sum of the regular 
and singular contributions with the singular contributions represented by the power- 
like distributions represented in Table 2. Finally, the piece of the cross-section that 
initially possesses the non-resonant structure, we consider in the conventional fashion. 

Now let us consider the interference contributions. In the case of the interference 
between the contributions with the single-resonant and the double-resonant structure, 
in the cross-section the non-regular factors (xj ± i0)~ fc appear instead of one of the 
BW factors (k is a positive integer). This factor is reduced via Sokhotsky formula 
to the sum of PV(xj)~ h and 5k-i(xj). Thus we obtain contributions that are similar 
to those which have been already discussed in the case with the double-resonant 
structure. The other cases of the interference are considered by similar fashion. 

4 The observable cross-section 

Having calculated the hard-scattering cross-section a, we proceed to the calculation 
of the convolution integral (2.1). In what follows we concentrate on the singular con- 
tributions only. (Recall that the regular contributions can be numerically calculated 
after the complete definition of the weights.) In the dimensionless variables integral 
(2.1) takes the form 

M 2 f 

a{s) = mj^i J dx ' (4 + x>) 9{x ' ~ Xmin) 0(rr '' x) ° {x>) ■ (4 - 1} 



16 



Here ^ = 2M (l+x/4), Vs 7 = 2M (l+x'/4), = 2M (1 +x min /4). The 0(x', x) 

is the flux function, 

0(x',x) = /3 e (4+x) 2(1 ^ e) (8+x+x')^" 1 (x-x')f- 1 



1 + 



4 + x ,x 2 



4 + x 



6(x - x) . (4.2) 



In view of a(x') ~ ^(aj' + ax +a 2 ) and ai + a 2 < 4, the integration area in (4.1) is 
restricted so that x' > —4. Consequently only the factors (x — x')^^ 1 and #(:r — x') 
are unsmooth in the integration area. 4 Owing to j3 e > both they are integrable. 

Further, we consider separately two cases, the first one when the value of external 
variable x is outside of the physical thresholds x = and x = —ai, and the second one 
when x falls precisely on one of the thresholds. In the former case the singularities 
of (x — x') 13 ^" 1 and 9{x — x') do not intersect with the singularities of a(x'). So 
all singularities are isolated and integral (4.1) can be numerically calculated with 
the aid of formula (3.13). In the case when x coincides with one of the thresholds, 
the singularities intersect. However, owing to the relations (±rr)+ — and 

(±rr)+ (±x)+ = (±:r)+ 1+A2 valid in the framework of the analytical regularization, we 
get either zero contribution or a single power-like distribution instead of the product 
of singular distributions. So integral (4.1) can again be numerically calculated by 
means of formula (3.13). 

Unfortunately, in the vicinity of the thresholds the behavior of the convolution 
integral is greatly unsmooth by the above procedure, and therefore the above de- 
scription cannot be considered satisfactory. Namely the integral (4.1) is unboundedly 
increasing with s approaching the thresholds (although remaining finite precisely on 
the thresholds). For definiteness let us consider the case x — » 0, which means s — > 4M 2 
in dimensional variables. In this case the dangerous singularities in a(x') are (±x')+, 
where 7 is half-integer, because with x — > the position of singularity of (x — x')+> 
where v is either f3 e —l or 0, is approaching that of (±x') + . At 7 + // < —1 this implies 
the increasing of the integral as Irrl^ 7 " 1 ^. The similar behavior takes place in the 
vicinities of the thresholds x = —ai, or s = (^s imin + Mj) 2 in dimensional variables. 

In effect the mentioned behavior means a rapid divergence of the asymptotic 
series when s is approaching the thresholds, and thus the expansion near thresholds 
becomes senseless. (Actually this behavior have already been detected in [28] where 
also it has been revealed that the area of bad behavior of the cross-section is rather 
small.) Nevertheless, the situation may be corrected by means of doing an additional 
expansion of the cross-section a(s) near thresholds, namely the expansion in powers 
of the distance between s and the appropriate threshold. For definiteness let us 
consider the physically important case when s is close to but a little bit lower of the 
threshold s = 4M 2 . In the dimensionless variables this means x < 0. In this case in 



4 In the strict sense, the factor 9(x' — a; min ) in (4.1) is unsmooth, too. However if a; min is the 
minimal kinematically- allowed energy, this factor already is present in a(x') and therefore can be 
omitted in formula (4.1). In the case x min > — (ai + a 2 ), we suppose that ir min does not fall on the 
physical thresholds. 
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a(x') the singularity (— x')\ is crucial. After changing the variable x' — ■> rr — £, the 
corresponding contribution to <r(s) takes the form 



Here ■?/>(£, x) is a regular contribution. Taking Taylor of (£— x)+ in powers of a; (there 
is no necessity to expand also the weight) gives the series 



instead of (£—#)+. Since the weight is regular in the vicinity of £ = 0, the integral 
of each term of this series can be numerically calculated. Truncating the infinite 
series, we obtain an asymptotic estimate of a at small x. Notice that this procedure 
provides a continuity of the result at deviating s from threshold. The number of terms 
in the truncated series should be chosen depending on the value of the deviation of 
s and of the current order of the expansion in a. Ultimately this number may be 
determined by the joining of the solutions outside and inside the interval where the 
second expansion of a(s) is made. 

So, outside of the physical thresholds the asymptotic expansion of the cross-section 
can be numerically calculated with the aid of formula (3.13). Near thresholds one 
should do the two-fold expansion, one in powers of a and the other in powers of the 
distance between s and the appropriate threshold. 



In this section we discuss the characteristic features of the application of the MPT 
within the NNLO approximation. For methodological reasons we consider the ex- 
pansion in the coupling constant in a sequential manner. By the first we consider 
the expansion of the BW factors, which by our definition let constitute the MPT ex- 
pansion in the narrow sense. From this point it is convenient to trace more carefully 
the field of subsequent calculations. Then we consider the expansion of the rest of 
contributions to the cross-section. 

First of all we consider contributions that possess the double-resonant structure. 
In this case the LO in the narrow sense, in accordance with (2.10), is determined by 
the substitution of the ^-functions (without the derivatives) for the BW factors in 
formula (3.2). Correspondingly, the $ must be determined by the on-shell amplitude 
multiplied by the branching factors, both calculated as a whole with the two-loop 
precision. The soft massless-particles contributions are to be taken into consideration 
by conventional fashion. In particular, at calculating the amplitude the Coulomb 
singularities are to be extracted in favor of the factor in the cross-section. However, 
contrary to the wide belief, this factor can be taken into consideration not in the form 
(2.13) only. In fact, it can be taken into consideration in the expanded form, as well, 




(4.3) 




(4.4) 



5 Calculations up to NNLO 
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though with certain modifications in subsequent calculations. Up to the NNLO three 
terms of the expansion are needed only, so we put 




2 



(5.1) 



The subsequent calculations can be carried out with replacing @ by (3 in the above 
formula, where (3 ~ \Jx — x\— x 2 is the off-shell velocity in the c.m.f. Since the process 
in the end is considered on-shell, this replacement should not affect the ultimate result. 
However at the intermediate stage of calculations it implies an effective modification 
of the kinematic factor. At the calculation of the hard-scattering cross-section, the 
second and the third term of (5.1) manifest themselves through the basic integrals 



A°! and A~\' 2 . From (B17) we have A^(x) = 6{x) and A~l' 2 {x) = x+'\ The 



convolution integral of such contributions without problems can be calculated. 

An alternative mode of taking into consideration Coulomb singularities has been 
discussed in Sect. 2. It implies a regularization of the inverse powers of (3$ in (5.1) by 
means of the special resummation. As applied to the second term in (5.1) this means 
a substitution 



The substitution of this kind is inevitable at the calculation in the NLO of the MPT 
in the narrow sense and in the higher orders, because in this case not only the on-shell 
contributions of the Coulomb factor must be taken into account but its derivatives, 
too. Fortunately, in the NLO in the narrow sense, but remaining within the NNLO in 
the general sense, it is sufficient to take into consideration only the one-photon contri- 
butions to the Coulomb factor, which is given by formula (2.15). Simultaneously, the 
weight $ must be taken into consideration off-shell and within the one-loop precision. 
The soft massless-particles contributions that are not involved in the Coulomb factor 
can be taken into consideration by means of technique discussed in [25]. 

In the NNLO in the narrow sense the weight $ should be taken off-shell in the 
Born approximation. The contributions of soft massless particles are absent in this 
approximation. Nevertheless, as long as the Coulomb factor makes considerable con- 
tributions in the threshold region, its contribution in the advance can be taken into 
account via formula (2.15). 

Further, we discuss the contributions that have single-resonant structure. At the 
amplitude level they arise either due to the process of inclusive single production of 
one of the unstable particles, or due to non-factorizable corrections in the mode of 
pair production. To the cross-section such processes can contribute in the pure form 
without the interference, or in the interference with the process of pair-production. 
Within the NNLO in the general sense the mentioned contributions can appear in 
the LO and NLO of the MPT in the narrow sense, and with the weight calculated 
within the NLO of the conventional PT and in the Born approximation, respectively. 
The unconventional counting of the orders is explained by an occurrence of an addi- 
tional power in a when one of the unstable-particle propagators squared disappears 
in the cross-section; cf. (2.5) and (2.7), and see discussion in [26, 27]. For the same 



2(3 



r.h.s. of (2.15). 



(5.2) 
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reason the interference between the non-resonant process and the single- or double- 
resonant process contributes to the NNLO in the general sense when the processes 
are calculated in the Born approximation. 

After the orders of the calculations are determined and the corresponding weights 
are calculated, the further calculations are to be made in accordance with the scheme 
of Sects. 3 and 4. Here we notice only that in the NLO in the general sense the 
following basic integrals are r equired: l\ /2 (x), A 1 /?, , B\ /2 , B 1 / 2 , b\ /2 . In the 
NNLO the additionally required basic integrals are j]^ 2 (x), J^ 2 (x), A 1 ^ 2 , A 1 ^ 2 , A]/ 2 , 
A 1 / 2 , C^ 2 , C\ 2 , C^ 2 , '■ The explicit expressions for them follow immediately 
from the appropriate formulas of Appendix B. The distribution sense is strongly 
required of the basic integrals A l J 2 \ A 1 / 2 , and C^ 2 2 . 

6 Discussion and conclusion 

We have constructed an asymptotic expansion in powers of the coupling constant a 
of the cross-section for pair production and decay of fundamental unstable particles. 
The algorithm of the calculation of the coefficient functions of the expansion is elab- 
orated in an arbitrary order. The outcomes are presented in the form of absolutely 
convergent integrals with precisely definite structure. Certain components of the inte- 
grals correspond to singular contributions in the amplitude and such components are 
process- independent. The MPT is intended for handling these singular contributions 
at the level of the cross-section, while the regular contributions are to be determined 
by means of conventional PT. The key moment of the MPT is the interpretation of 
the singular contributions in the distributions sense, so that they become integrable. 
In doing so, the asymptotic property of the expansion is strictly maintained. The 
exception is made for the Coulomb singularities, which are considered extracted into 
a factor and represented in the regularized form in the sense of conventional func- 
tions. However the regularization implies no insertion of extrinsic parameter but 
special resummation of intrinsic contributions [6]- [9] So the asymptotic property of 
the expansion in full measure is maintained. 

In effect, the difference between the conventional PT and the MPT appears at 
calculating the cross-sections and becomes apparent at taking into consideration the 
resonant contributions of the unstable particles. At comparing the MPT with the 
other approaches intended for the description of the production and decay of unsta- 
ble particles, an essential moment distinguishing the MPT is the complete expansion 
in powers of the coupling constant a of the resonant contributions of the unstable par- 
ticles. In view of this property one can expect the gauge cancellations at calculating 
the cross-section in the framework of the MPT approach. 

The MPT method in the equal measure is applicable for the calculation of the 
total cross-sections, the angular distributions, and spin correlations. This follows 
immediately from the fact that the angular and spin variables are not involved in 
the BW factors. The calculation of the invariant-mass distributions is workable in 
the MPT, too, but in conjunction with some additional operations. Namely since 
the BW factors in the MPT do not have the sense of conventional functions, an 
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additional integration of the BW factors with some weight is necessary. In any case the 
integration should be made in view of the jets and the contributions of unregistered 
photons in the final state. Unfortunately, the convergence property of the MPT 
becomes worse at the proceeding to the invariant-mass distributions. For this reason 
we do not suggest the use of the MPT in this case, at least in the considered above 
mode. For similar reason we do not suggest the use of this mode of the MPT for the 
description of a generic 2^2 process with a single resonance. Instead in these cases 
it is better to use a modification of the MPT approach discussed in [26, 27], which is 
based on the secondary Dyson resummation in the framework of MPT. 

The question of the great importance in the case of pair production of unstable 
particles is the near-threshold behavior of the cross-section. Fortunately, this behav- 
ior can be calculated in the MPT framework, as well, though through the double 
expansion of the cross-section. Namely apart from the ordinary MPT expansion in 
powers of a, one should do expansion in powers of the distance between s and the 
appropriate threshold. The necessary number of terms that must be retained in 
the latter expansion depends on the value of the distance and on the current order 
of the expansion in a. Ultimately it is determined by the claimed precision of the 
description. 

The convergence properties of the MPT in the case of pair production of unstable 
particles have been preliminaryly studied in the framework of model calculations, in 
the above-threshold region within the NLO [28]. (The modelling has concerned the 
weights, but not the structure of the MPT itself.) The similar calculations within 
the NNLO and with taking into consideration the Coulomb singularities have been 
recently made, as well [48]. In both cases a satisfactory behavior of outcomes has 
been detected with essential improving of the description at the proceeding to the 
NNLO. Numerical analysis in the threshold region will be carried out elsewhere. The 
calculations of the completely realistic processes will follow. 

In conclusion, the basic result of the given work is that the very existence is proved 
of an asymptotic expansion in powers of the coupling constant of the cross-section 
for pair production and decay of fundamental unstable particles. Furthermore, a 
scheme of the calculation of the expansion in an arbitrary order has been explicitly 
constructed and its specific ingredients have been calculated. On the whole, the MPT 
method turns out to be a real candidate for carrying out high-precision calculations 
necessary for the realization of programs at colliders of next generation. 
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Appendix A 

Here we derive formula (3.4) for the change of variables in the 5-function with 
derivatives. More precisely, we derive a formula for the distribution 5 n (p(x)) resolved 
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with respect to x variable with 

p( x ) = x + bx 2 , b^O. (Al) 

At once we note that since the distribution S n (p) is concentrated at p = (i.e. gives 
zero contributions outside a neighborhood of p = 0), and since at the inverse mapping 
p — > x the point p = turns into two points x = and x = —1/6, the sought-for 
formula generally has the form 

n n 
fc=0 fc=0 

Here Zfc and Cfc are coefficients. However, our actual interest is with the formula 
considered in a neighborhood of x = 0. For this reason we reduce the problem to the 
calculation of the coefficients z k only. 

First we note that in the case of small enough neighborhoods of x = and 
p — 0, formula (Al) defines the one-to-one mapping p <-> rr, with the inverse mapping 
determined by formula 

^v/TT^-l (A3) 



Further we note that by virtue of 

dxx n 5 m (x) = 5 nm (A4) 

and in view of (A2), the coefficients z k are determined by the formula 

z k = J dxx h 5 n (p(x)) , (A5) 

where < k < n and the integration is carried out over the above mentioned neigh- 
borhood of the point x — 0. Making the change of variable x — ► p, we get 

An elementary calculation based on the formula for binomial differentiation yields 

-ffisO(-^W"^) r (^) r ( i+ -^)- (A7) 

In this formula at odd r the sine under the sum gives zero contributions. So we can 
set r = 2m, m = 0, . . . , [k/2], where [k/2] is the integer part of k/2. After further 
manipulations with the gamma-functions, we obtain 

_ ^)k+n b n k] [m (_yn ( 2w -2m)! 

Zk ~ (2b) k n\ ^ m\ [k - 2m)\{n - m)\ ' ( ^ 
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The calculation of the sum yields [49] 



z k = (^^ {-bY~ k . (A9) 

Now we substitute (A9) into formula (A2) considered without the second sum. 
Making the change of variables k — > n — k, we obtain 

S n (p(x)) = £ ( n+ n k ) (-*)* ^-k{x) . (AlO) 

Let us remember that the contributions located at x = —1/b are omitted in this 
formula. 

A generalization to the case of presence of a common scale factor in formula (Al) 
is trivial and is practically implemented in (3.4). 



Appendix B 

In this Appendix we analytically calculate the single-subscript basic integrals 1^, 
and the double-subscript basic integrals A£ , B^ in2 , C* in2 , considered with posi- 
tive parameters Oj and non-integer A. In the single-subscript integrals that depend on 
one of the parameters Oj, we omit the lower index "i" for simplicity of the notation. 
In the case of the double-subscript integrals we introduce a = a\ + a 2 and n = n\ +n 2 . 



Calculation of \\ 

Let us discard the superfluous variable Xj in formula (3.16) and omit the index 
"i" in the parameter a«. Then we come to an equivalent definition, 



I x n (x) = J d£0(£ + a)(s-Oi<5»-i(O- (Bl) 

Here n is a positive integer, a and A are continuous parameters, a > and A is 
non-integer. Recall that if A is large enough, the (x — means 6(x — £) (x — £) A . 

At first stage we consider A as large enough. By making the change of variable 
£ — > x — £ and resolving the ^-functions, we come to the equivalent formula for the 
integral written in finite limits: 

= ^* (B_1) « " *) • (^2) 

Notice that at rr < integral (B2) vanishes owing to the strict positiveness of the 
argument in the 5-function. At x > the calculation by parts yields 

ll{x) = e ^±^{-T-H{x){x^ n - l) . (B3) 
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Further we can omit factor 6(x + a), because at a > the condition x + a > is 
automatically satisfied in the presence of 9(x). After calculation of the derivatives, 
we obtain 

It is readily seen that at non-integer A this expression is well defined for any x 
except x = 0. At x = the above expression, generally, is ill defined. However 
eventually we do not need a value of the function I x (x) at x — 0, but a rule of its 
integration in a neighborhood of rr = with smooth weight. In the framework of 
analytical regularization the 6(x) x 1+x ~ n means the distribution x+~ x ~ n . So in the 
case of analytical regularization we have 

With transiting to another regularization the result may change. To make this 
more clear, let us consider an alternative derivation of formula (B5) by basing on the 
direct calculation of the integral in (Bl). In this way at first we formally obtain 

1 rT -1 

= ^ H + *M* ~ (* " 01 1 ^ ■ (B6) 

Since with a > the differentiation of #(£ + a) gives zero contribution at £ = 0, we 
omit this ^-function. The differentiation of other factors gives 

= feS * w + S § (V) ^) l "" r_1) ■ < B7 > 

In this formula the sum by definition equals zero if n — 1. If n > 1, the contributions 
under the sum generally are ill-defined and so they need a regularization. Generally, 
the different regularizations lead to different outcomes. For instance, under the cut-off 
regularization the sum takes the form of Ylr=i ^ r ^ r \ x )i where C r are non- vanishing 
coefficients. Under the analytical regularization all coefficients C r vanish. Really, in 
the framework of the analytical regularization formula (B7) must be considered at 
first at A > n — 2, when all the terms in the r.h.s. are integrable. In this case all 
terms under the sum become zero in view of the presence of S^ix). So only the 
first term in (B7) survives. Having calculated the derivatives in this term, we come 
to formula (B5) again. 

Calculation of J A 

Proceeding on analogy to the previous case, we rewrite definition (3.17) in the 
form 

3 x n (x) = J d£ Oft + a)(x- Oi PV± . (B8) 
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Then, making the change of variable £ — > rr — £ and resolving the ^-functions, we get 

= e(x + a) (-)» £ +a d£ e PV j^r ■ ( B9 ) 

Further let us consider an auxiliary integral with b>0, 

Jn\b,x) = ^ d ££ A P^-^. (BIO) 

At A > n — 1 and x ^ 0, x ^ b this integral can be calculated by means of reducing 
the improper integral to the derivatives of the absolutely convergent integral: 

C b 1 1 d n f b 

^ - - ^ W„ d? « A i« -I- (bid 

The integral in the r.h.s. can be explicitly calculated. At b > x and non-integer 
A > —1, we have 

' b b l+x , , i .,/_. , . . or 

b, 



r b b 1+x r 1 / 

J d££Mn|£-x| = — In(6-x)-— J^l,-A-1;-A; 



(B12) 



+ F(1 + A) ^ A) [-cos(n\)6(x)x 1+x + 6(-x) (-x) 1+A ] . 
1 + A 

Hereinafter F means the hypergeometric function 2 i*i [50]. Substituting (B12) into 
(Bll) and carrying out the differentiation, we get 

Jn(b,x) = — — F(n,n-\-l;n-\;-) (B13) 

1 + A — n \ o/ 

1 (raj 

Comparing (B9) with (BIO) and (B13), we come to the formula: 

J x (x) = (-)" {x + a) } +X - F (n,n-\-l;n-\;—) (B14) 



1 + A — n \ x + 

+ 9(x + a) r(1 + A) ^ ) ~ A ~ 1) [- cos(ttA) x\ +x - n + (-)» (-x) 1 ^] . 

We see that J x (x) has, in general, two singular points, the x = and x = —a. 
At x = the singularity is controlled by the second term in (B14). To find the 
behavior at x — > —a, we take advantage of the analytic-continuation formula for the 
hypergeometric function with inverse argument [50]. In this way at x < we obtain 

x/ \ (x + a) 1 + x / , , x + a 
J n(^) = / Wl xx i 71 ^, 1 + A; 2 + A; 



(-x)»(l + A) V a; 

(B15) 

This formula is equivalent to (B14) at — a < x < 0. With 1 + A > the above 
expression is regular at x + a — > 0. 
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Calculation of A* iH2 



Let us consider integral (3.10) as the iterated one, and let us at first calculate 
integral in x 2 . With the aid of (Bl) and (B5) we obtain 

<n» = r(njr(2 + ^) I dXl 6{Xl + ai) {X " Xl) + +A " n2 6nMXl) ■ (B16) 
The similar calculation of the remaining integral yields 

A A (r\- (-)" r(l + A) 2+A-n mi ? A 

A n in2 W - r(ni) r(n2) r(3 + A _ n) * + • (B17) 

Here n = n\ + n 2 . The result is symmetric and does not depend on the (positive) 
parameters a^. 



Calculation of B ^ 



As in the previous case, we consider integral (3.11) as the iterated one and at first 
calculate integral da; 2 . With the aid of (Bl) and (B5) we get 

= r^SK / dXl Hxi+ai) <* - PV '^ • (B18) 

Then, with the aid of (B8) and (B14) we obtain 



(X 
n 1 ,n — X — 2;n — X — l; 
x + ai 

~ 0{x + fll ) r(1+ ^ )r , ( r 7^" 2) [cos(ttA) x^ + (-)» • 
i (niji (n 2 J 

The result is independent from the order of the calculation of integrals, which may be 
verified by direct calculations. Notice also that (B19) is independent from the param- 
eter a 2 . The conjugate integral B ^ iri2 is derived by the substitution {ni,n 2 ,ai} — > 
{n 2 ,ni,a 2 }. 

In general, there are two singular points in (x), the x = and x + a\ = 0. 
At x — > the singular behavior is completely controlled by the second term in (B19). 
The behavior at x + a\ — > can be found on the basis of the representation of the 
hypergeometric function F in the form with inverse argument [50]. At x < we 
obtain: 

ni M , Nn _! T(l + A) (x + fl!)^™ 2 ^( n . , x + ai' 



r(n 2 ) T(3 + A - n 2 ) (-x)" 1 \ x 

(B20) 

At 2 + A — n 2 > this expression is regular as rr — > —a\. At 2 + A — n 2 < it contains 
singularity (x + ai)+ +A ~™ 2 and also the associated set of singularities (x + a 1 )^~ x ~ n2+r , 
where r is a positive integer such that 2 + A — n 2 + r<0. 
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Calculation of C* 

niri2 



The calculation of C„ ira2 is a more complicated task. We begin with the ob- 
servation that at large enough A the singularities in the integrand in (3.12) do not 
intersect. So the integral is well defined as the iterated one and theoretically can be 
calculated on the basis of the above results. However in view of complexity of the 
intermediate expressions, the analytical calculation in this way unlikely is efficient. 
An efficient method from the very beginning should allow for symmetrical proper- 
ties of the integral, which are evident when considering the integral as the double 
one. (Notice that a double integral of the PV-poles and of the 5-functions may be 
defined through the representation of these singular functional as improper limits of 
the regular functionals, see for instance [30] and [53].) 

So let us consider integral (3.12) as the double one. For the reasons that be- 
come clear later, we begin calculations with adding and subtracting a certain term 
proportional to A£ . Namely we rewrite (3.12) in the form: 

Cn in2 (x) = {-) n+l Tx 2 A^ in2 + //dm dx 2 6( Xl + ai) Q[x 2 + a 2 ) {x-x l -x 2 ) X + 



x i pv -^r + (-r^K-iW^-iM \ . (B2i > 

x 2 



Note that the integration in (B21) goes over a simplex. So for providing symmetric 
integration, we make a shift Ob >i ^ Ob 2, Q>% 

and then a transition to the cone variables 
x\ + x 2 = £, x\—x 2 = 2rj. After that, formally representing the double integral again 
as the iterated one, we obtain 

x+a 

C n \ n2 (x) = (-r + V< n2 + 9(x + a) (-rj dt(x + a- £) A (B22) 

o 

x I dr]\ PV- ^ — PV- 1 



(rj + C/2 - ai )"i ( V -i/2 + a 2 ) 



+ (-r+V 5^(7] + i/2 - ai ) 5^(7] - i/2 + a 2 )}. 



Unfortunately, in this formula the expression in the curly brackets, generally, is ill- 
defined. Really, with x > the variable £ may become equal to a = a\ + a 2 with the 
consequence that there appear the intersecting singularities. However this behavior, 
actually, is an artefact of the transition to cone variables and the problem may be 
removed on the basis of the following relation [51] (which, in turn, is a generalization 
of the relation obtained in a simpler case [52]): 

pv j^r^i PV X^r + ( " rV 6ni - l{z " Zl)6n2 - l{z " Z2) (B23) 

711-1 1 

= E rr 1 ) (-ypyjT—^py 
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«2 - 1 
r=0 



^2 - 2l)" 1+r 



The meaning of this relation is as follows. Let us assume that some integral with a 
weight of the expression in the l.h.s. is to be calculated at first in z and then in z\ or 
z-i- Then the calculation of the integral in this sequence should be carried out with 
the substitution of the expression in the r.h.s in this formula. In our case, we need 
a projection of the relation (B23) obtained by substitutions z = rj, z\ = — £/2 + a±, 
z 2 = £/2 — a,2 ■ After the substitutions, the l.h.s. in (B23) becomes the same expression 
as the in the curly brackets in (B22). So the initial integral can be converted to the 
form of (B22) with substituting the following expression for the expression in the 
curly brackets: 

m-i 1 1 

V^na+r-n (_)rpy py " ( B24 ) 



r=0 
n 2 — 1 



a _£)n 2 +r * ' ( ?7+ £/2-ai) ni - r 

1 



(?7-e/2+a 2 )^- r ' 

A consideration of completely solvable example [51] independently shows that the 
substitution of (B24) solves the problem of the ghost singularities emerging at the 
transition to the cone variables. 

So in the integral di] in (B22) we gain the contribution of a sum of isolated PV- 
poles instead of intersecting singularities. The integral of the poles, we calculate with 
the aid of the formulas 

,6 1 

drjPV- = ln|fc-£| -ln|a-£| , (B25) 



I 

J a 



rj 



b 

I 



dr] PV- 



(v - n 



n- 



PV 



1 



(b - o- 1 



PV- 



(a - 



(B26) 



After the calculation of the integral d.77, we make the change of variables £ — > x + a — £ 
in the remaining integral d£ and simultaneously the change r — > n« — r — 1 in (B24). 
This yields (a = a\ + 02) 

(A 1 „\ 1 \n+l 2 4 A 



(B27) 



+ 6{x+a)- 



\n-l 



r(n-i) 



r(ni) r(n 2 ) 



| d£ £ A (ln 


£ — a; — ai 


+ In 


£-x-a 2 






a 2 



PV- 



n-l 



+ 0( x + a)(-) 



n-l 



m — 1 

E 

. r=l 



_ 1 -1 x+a 

("„7-"i 2 ) +E C»7-"i 2 ) T" /d« A PV- 



ni — 1 



r 



(£_ x )n-r-l 
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719-1 X + CL 

1 r 1 1 

In this formula if the upper bound in any sum is less than the lower bound, then by 
definition the sum is zero. 

Let us discuss the obtained outcome. First we note that, actually, we have cal- 
culated the integral in the third term in (B27) while calculating the basic integral 
J*. The integrals in the two last terms can be easily calculated, as the singularities 
of the PK-poles do not coincide at a, ^ 0. So the only problem integral remains in 
the second term, which contains logarithm contributions. To calculate it, we subtract 
and add to each logarithm the n — 1 of the first terms of its expansion in powers of 
£ — x, so that to represent the logarithm in the form 

♦SK 5 ?)'}-^)'- 

Here the contribution in the curly brackets is of order 0((£ — x) 71 ' 1 ) as £ — x — > 0. 
Therefore the appropriate contribution to the integral is regular in the vicinity of 
£ = x. The contribution of the last term in (B28) is reduced to the sum of integrals 
considered at the calculation of J*. 

After carrying out all above mentioned calculations, we come to the following 
result: 

c n \ n2 (x) = (-r + v<„ 2 (*) + ^^rin^ I d ^( x+a -oe + xn^x), 

(B29) 

+gH^r)' + (™)}' (B30) 

Since Xn((,,x) by the construction is a regular function, integral (B29) is absolutely 
convergent at A > — 1. In this case the singularity of C^ in2 (x) is contained only in 
the A^ iri2 (x) which is present in (B29). 



In 



£ - X - di 


= jln 


£ — X — di 

















n— 1 



In 



£ — x — di 



di 
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